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The explicit expression for the Luttinger number N L corresponding to the SU(iV) model of Dave, 
Phillips and Kane reveals that N L is equal to the number of particles n only when n = N/2 (for iV 
even and assuming that the unit-step function 0(x) is equal to ^ at x = 0) and n = N, signalling 
failure of the Luttinger theorem for all other values of n in the interval {1, 2, . . . , TV}. In this 
Comment, we first present general arguments showing that the absence of the Luttinger theorem 
for the SU(A) model under consideration is rooted in the non-uniqueness of the ground state of 
this model for < n < N, the validity of the Luttinger theorem for n — N/2, when TV even, 
being accidental, a consequence of particle-hole symmetry. Consequently, by supplementing the 
Hamiltonian of the SU(7V) model with a perturbation Hamiltonian that removes the ground-state 
degeneracy, the Luttinger theorem is to apply for the resulting model in the limit of the coupling 
constant A of this perturbation approaching zero, where the limit A — > is clearly to be taken 
subsequent to taking the zero-temperature limit of the thermal single-particle Green function in the 
expression for N L . We explicitly establish the validity of this statement for the case of N = 4. The 
details of the relevant calculations being distinctly transparent, one can readily convince oneself 
that our observation is valid for arbitrary N. It follows that the issues raised by Dave, Phillips 
and Kane, such as non-existence of the Luttinger-Ward functional and "breakdown of the elemental 
particle picture in strongly correlated electron matter", are all inessential to the observed failure 
of the Luttinger theorem. As regards the singularity of the self-energy £(u)) on the real w-axis, 
observed by Dave, Phillips and Kane, we demonstrate that this also is a direct consequence of the 
non-uniqueness of the ground state of the SU(iV) model for < n < N. In the light of the above 
observations, we are in a position to state that to this date no case has come to light indicative of 
the failure of the Luttinger theorem under the conditions for which it has been deduced. 

PACS numbers: 71.10.-w, 71.10.Pm, 71.27.+a 
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I. INTRODUCTION 

In a recent publication, Dave, Phillips and Kane (here- 
after DPK) [1] have considered the many-particle system 
described by the Hamiltonian 

n = ^(h)\ (i) 



where U > is the interaction-energy parameter, and 

N 

n=^2 "«i where n a = c^c a , (2) 

<x=l 

in which {c a } are canonical annihilation fermion opera- 
tors, and {etc} their Hermitian conjugates. 

The thermal single-particle Green function ^ aa ' (uj) 
pertaining to the SU(iV) model under discussion proves 
to be diagonal with respect to the indices a, a' [1] . Con- 
sequently, with G aa '(a;) denoting the zero-temperature 
limit (corresponding to (3 = 1/T — >• oo - throughout 
this Comment k B = 1) of the ^^(w) specific to the 
zero-temperature limit of the chemical potential fip cor- 
responding to the mean value of particles n in the grand- 
canonical ensemble (see Sec. II Bl however), the Lut- 
tinger number (LN) N L [2, §2] is denned as follows [3, 
Eq. (95)]: 

N 

N L = ]T©(G aa (0)). (3) 

a=l 

According to the Luttinger theorem (LT) [2-4] , one must 
have 

A^l = n. (4) 

For the SU(JV) model, DPK [1] have however deduced 
that 

A L = N0(2n — N). (5) 
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Assuming that 0(0) = \ (see however Rcf. [2, §2.4]), one 
observes that unless n = N/2 (for N even) or n = N, 
the LT is clearly violated for the SU(AT) model under 
consideration [1]. In this Comment we clarify the reason 
underlying this failure. Briefly, by the non-uniqueness of 
the ground state (GS) of the SU(iV) model for < n < 
N, this model is for < n < N a priori excluded from 
the set of models to which the LT is applicable. That 
despite the apparent non-uniqueness of the GS of this 
model for < n < N the LT proves to be valid for N 
even and n = N/2, is a direct consequence of the explicit 
particle-hole symmetry of the problem at hand when in 
this case the chemical potential \i is identified with the 
zero-temperature limit of 

II. ANALYSIS 
A. General considerations 

Let \\Po.i) = |0,0, . . . ,0) (N zeros) denote the normal- 
ized vacuum state of the SU(iV) model under consider- 
ation. For the normalized n-particle eigenstate \& n ,s) of 
rl, with n > 1, one has: 

\#n,.)=n? = i&l U) \*o,i), SG{1,2,...,(J]}, (6) 

where the integer- valued mapping 

v s :{l,...,n}^{l,...,N}, (7) 

for a given s, coincides with an TV-permutation Vn.s' 
of the ordered set {1, 2, . . . , N}. Here, the two per- 
mutations Vn.s' and Vn,s" are to be identified when 
the ordered set {'Piv.s'l, . . . , Vn,s'Ti} is an n-permutation 
of the ordered set {Vn !S "1, ■ ■ ■ ,Vn, s "Ti}. Equivalcntly 
as regards the ordered sets {VN,s'{n + 1), . . . ,Vn,s'N} 
and {VN,s"{n + 1)> • • • > Vn,s"N}. Hence, one is left 
with N\/(n\(N — n)\) = ( ) distinct permutations, a 
fact reflected in the specific set over which the inte- 
ger s in Eq. (6) varies. Note that because of the anti- 
commutation relation [c ?c«']- = 0, Voc, <x', identifying 
v s with either Vn,s' or Vn,s", with {Vn,s'1, ■ ■ ■ ,VN,s'n} 
and {'Pjv>"1, ■ ■ ■ , 'PN,s' , n} differing by an n-permutation, 
the corresponding \& n ,s) differ by at most a minus sign. 

For later reference, the state |^ ra , s ), Eq. (6), can be rep- 
resented as follows (compare with the l^o.i) introduced 
above) : 

\&n,s) = |n S) i,n Si2 , ■ ■ .,n a ,jv), (8) 
where n SiCt £ {0, 1}, Va, and 

N / A A 

5>,« = *V,e{l, 2... .,( )}. (9) 

From 

H\& n , s ) = ^n 2 \& n , s ), Vs, (10) 



it is observed that for n particles the energy level Un 2 /2 
of the SU(iV) model is ( )-fold degenerate. For later ref- 
erence, we point out that the GS energy of the Hubbard 
Hamiltonian [5, 6] for n spin-^ particles defined on a lat- 
tice comprised of N sites is in the atomic limit [6, §5] 
2 n (^)-fold degenerate for n < N [7, §4.4.2], where the 
factor 2™ corresponds to the spin degeneracy; at half- 
filling, corresponding to n — N , one is thus left only with 
spin degeneracy. For definitencss, with {Tj^} denoting 
the hopping integrals in the Hubbard Hamiltonian, the 
atomic limit is defined by the condition Tj j = for all 
i 7^ j [6]. The term T,^ = To is arbitrary and is usually 
absorbed in the chemical potential fj, (cf. Eq. (17)), how- 
ever it has been explicitly taken into account in Rcf. [6]. 

Two main problems arise as a result of the degeneracy 
of the GS energy of the SU(iV) model for n ^ 0, N. 

First, in the cases where the GS is not unique, for fi G 
>/•*«)' Eq. (33) [2, Eq. (B.7)], the zero-temperature 
limit of the thermal single-particle Green function coin- 
cides with an average over the set of zero-temperature 
single-particle Green functions each member of which 
corresponds to one of the GSs in {Ji^s) \\s}. This av- 
erage function has no place in the formulation of the 
LT. In Ref. [2, Appendix C] we have at places indicated 
the consequence of degenerate GS energies (see for in- 
stance Eq. (C.22) and the remark following Eq. (C.23) 
in Rcf. [2]; the expression in Eq. (C.22) is to be viewed 
in the light of the expressions in Eqs. (C.12), (C.19) - 
(C.21)). Later, in Sec. IIB2, we explicitly show how 
the zero-temperature single-particle Green function ob- 
tained through the thermal averaging over the manifold 
of the single-particle Green functions corresponding to 
the n-particle GSs {|#n, s )||s} in general yields an incor- 
rect (from the perspective of the LT) LN 7V L , Eq. (3). 

Second, the LT has been explicitly proved within the 
framework of the time-dependent many-body perturba- 
tion theory [8] (for a comprehensive discussion, see in 
particular Ref. [2]). In this framework, it is assumed 
that at zero temperature an adiabatic increase of the 
coupling constant of the interaction Hamiltonian, from 
zero at t = — oo to the full strength at t = 0, results 
in an adiabatic evolution of the unique (up to a triv- 
ial phase factor) n-particle non-interacting GS into the 
unique 7i-particle interacting GS. The latter state is fur- 
ther assumed adiabatically to evolve, on adiabatically 
decreasing the strength of interaction to zero, into an 
7i-particle state at x = oo that up to a possible trivial 
phase factor coincides with the non-interacting n-particle 
GS with which one has begun at x = — oo. According to 
the celebrated Gcll-Mann and Low theorem [9] , on which 
the time-dependent many-body perturbation theory is 
based, the first of the above-mentioned adiabatic pro- 
cesses in general generates an n-particle eigenstate of the 
full Hamiltonian at x = 0, not necessarily its n-particle 
GS [8, p. 61]. Failure of the state at x = to be the 
interacting GS, amounts to a failure of the many-body 
perturbation expansion [10] (see also in particular the 
remarks in §5.1.1 of Rcf. [2]). 
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For completeness, as demonstrated in Rcf. [2], for a 
class of systems, fully described in the latter reference, 
many-body perturbation theory - in the specific way it 
has been employed in the proof of the LT [4], does not 
break down. Briefly, expansion in terms of the skeleton 
self-energy diagrams [4] [2, §5.3.1] evaluated in terms of 
the exact interacting single-particle Green function safe- 
guards the relevant perturbation series against either di- 
verging or converging to false limits. 

The LT having originally been implicitly deduced for 
metallic GSs, on account of an earlier observation by 
Rosch [11], in Ref. [2, §6.1] (see also Ref. [12]) we have 
demonstrated that a possible failure of the LT for non- 
metallic GSs, when the chemical potential [x £ (/i~,/i+) 
differs from the zero-temperature limit of fip, is not due 
to a possible failure of the many-body perturbation the- 
ory, but due to the possibility of arriving at a false limit 
[13, §302-306] for N L in the process of effecting the zero- 
temperature limit {3 — ¥ oo for any fixed value of fi in- 
side (fj,~ , fj,n) different from (j,^; we have also demon- 
strated that the last-mentioned false limit is similarly 
avoided by identifying /i with \xp in the process of effect- 
ing the zero-temperature limit [14]. We note in passing 
that in Ref. [15] we have rigorously dealt with the half- 
filled GS of the one-dimensional t-t'-V model, which for V 
greater than a critical value V c (t,t') is a charge-density- 
wave state and insulating. For this insulating GS, the LT 
proves to be trivially satisfied for all \i inside the relevant 
gap (fi~,[i+) [15, §IV] (sec Sec. C 1 in appendix C). 



B. Explicit treatment of the SU(iV) model for N = 4 

1 . Removal of the GS degeneracy 

Below we explicitly consider the specific case of N = 
4, for which the relevant Fock space is spanned by the 
following set of J2 n =o (n) = ^ simultaneous eigenstates 
of H and h (cf. Eqs. (8) and (9)): 



|%l> = |0,0,0,0), (11) 



|#i,i) = |l, 0,0,0), |# x , a > = |0,1,0,0), 

|#i l3 ) = |0,0,1,0), ^,4) = |0,0,0,1), (12) 



*2,1> = 


11,1,0,0), 


1*2,2) = 


[1,0,1,0), 


#2,3) = 


[1,0, 0,1), 


1*2,4) = 


10,1,1,0), 


*2,5) = 


10,1,0,1), 


1*2,6) = 


|o, 0,1,1), 


*3,l) = 


[1,1, 1,0), 


1 ,7V \ 

1*3,2) = 


|1, 1,0,1), 


#3,3) = 


|i, 0,1,1), 


1*3,4) = 


10,1,1,1), 




I*m) = 


1,1,1,1) 





(13) 



(14) 



(15) 



In assigning an s to a particular state corresponding to a 
given n, where s £ {1, 2, . . . , (^)}, we have had in mind 
the following 'perturbed' Hamiltonian (cf. Eq. (2)): 



N 

ri\ = H + A ^ oth c 

tx=l 



(16) 



Two significant properties of this Hamiltonian are: firstly, 
that for all A £ R its eigenstates coincide with those of H 
and n, presented in Eq. (6), and, secondly, that its GS is 
unique for all A > 0. With reference to the Hamiltonian 
in Eq. (16), the energies of the states in Eqs. (11) - (15) 
are non-decreasing for increasing values of both n and s. 
Thus, for any n £ {1,2,3,4}, !?Vi,i) is the n-particle GS 
of U\ for all A > 0. 

Let {K s (n, A)} denote the eigenvalues corresponding to 
eigenstates {Itf'n.s)} of the thermodynamic Hamiltonian 



1C\=r{\- (lit, 



(17) 



where /i is the chemical potential. For the corresponding 
thermal single-particle Green function c S 0ia j ( w , A), one 
has (throughout this Comment, h— 1) [2, Eq. (C.19)]: 



N 



^ W ' A ' = m£ E e-^"< A >/« a ,( W ,A;n, S ), 

^ > n=0s£6(n) 

(18) 



where [2, Eq. (C.9)] 



N 



= £ E 



-pK s (n,\) 



(19) 



n=0 s€6(n) 

is the grand partition function, 

{!,..., (£)}, n<N, 
0, n> N, 

where denotes the empty set, and [2, Eq. (C.20)] 
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(20) 



foca>(u,\',n,s) 



E 



(A«<(n)) 5 y (Ai(n)) s >, 



E 



Jx(™+l)) S ,s'( A a'(™ + l))s', S 



uj - K s (n,\) + K s ,(n-l,\) ' ^£7^. w — K s '{n + 1, A) + K s (n, A) ' 



(21) 



in which 



(A a (n)) s 



I 

For the states {|^n,s)} given in Eqs. (11) - (15), one read- 



on- 



ly/ 



(22) 
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ily obtains (below, S a 



A«(l) 




A«(2) 



/5_oc,2 







^ 


$_a,3 





Sa.,1 





Oa.,4 








5<x.i 





<5<x,3 


Sa.,2 








<5«,4 





b~oc,2 


V 





Oa.,4 





A«(3) 



/ $0i,3 


5oc,2 





<5a.l 





\ 


0~a,4 





<5a,2 















#a,3 










V 








Soc,4 


<5«,3 


Sa.,2/ 



^a(4) = (5<x,4 #«,3 <5a,2 ^a,l)- 



(23) 
(24) 

(25) 



For the considerations of the present Comment, we need 
to calculate the function 



Gaa> {oj, A) = lim Sf«*/(w,A), A > 0, (26) 



corresponding to the n-particle GS of 7^ a, A > 0. 
From the above expressions, one readily deduces that 
^aa'i^J, A) is identically vanishing for cc ^ oc', V/3, A, sim- 
ilar to the case of A = considered in Ref. [1]. 

Since the n-particle GS of T-L\ is unique for all A > 0, 
equal to ^n,!}, following the asymptotic expression in 
Eq. (C.12) in Ref. [2] for fi e (^t~,/i+), one immediately 
obtains 



G aa (u;,A) = ^ 

see(n-i) 

+ E 

see(n+i) 



|(A«(n))i 



to - K^n, X) + K s (n - 1, X) 



V a (?l + l)) s 



u-K s (n + l,X)+K 1 (n,X) 



(27) 



Making use of the representation in Eq. (27) and the 
above expressions for the matrices {A a (n)||n}, one read- 
ily arrives at the following results: 



G« a (0,A)|„ =1 = 



J«.2 



H — U/2 — X fi — 3C//2 - 2A 

+ ±:i (28) 



H - 3U/2 - 3A ^ - 3C//2 - 4A ' 



G««(0,A)| 



J a., 2 



n = 4 H-7U/2-X H-7U/2-2X 

H-7U/2-ZX fj, - 7C//2 - 4A ' 
We introduce [2, Eqs. (B.5) and (B.6)] 

H~ = H{n) - H{n - 1) = (n - ijET", < n < N, 



(31) 



M+ = H(n + l)-H(n) 



(n+-)C7, n <N, 



(32) 



oc . 



n = N. 



where H (n) = Un 2 /2 is the degenerate eigenvalue of H 
corresponding to \& n ,a), s G {1, . . . , (T)}, Eq. (10). Note 
in passing that in contrast to Ref. [I], here we do not 
consider H(n + 1) as being a definite quantity for n = N 
(cf. Eq. (20)). It is observed that in the limit of A = + 
for [2, Eq. (B.7)] 



(33) 



the denominators of the first n terms on the right-hand 
sides of the above expressions for G aoc (0, A), with A = + , 
are positive. For TV = 4, one thus has the following exact 
result: 



N 

E 



©(G aa (0,0+)) =n, V M gfc,^), 



(34) 



implying the validity of the LT for all n S {1, 2, 3, 4}, not 
only for /i identified with the zcro-tcmpcraturc limit /ioo 
of the chemical potential fj.p corresponding to n particles 
(which in the case at hand is equal to Un [1] for n < N, 
and not less than Un for n — N), but also for all /i 
in the single-particle excitation gap (/z~,/U+) of the n- 
particle GS of the model under consideration, with A = 
+ . The equality in Eq. (34) is to be contrasted with that 
in Eq. (5). We have therefore rigorously demonstrated 
that in the specific case of N = 4 the failure of the LT 
as observed by DPK [1] is wholly attributable to the n- 
particle GS ofH\ = o, to be distinguished from %\=o+ » n °t 
being unique for < n < N . 



2. A consequence of the thermal averaging of the Green 
functions corresponding to different GSs 



Sa.,1 



G«a(0,A)| n=2 - fi _ m/2 _ x ■ /i _ 3C/ /2_2A 

<5«,3 StxA 



/X-5C//2-3A //-5Z7/2-4A 



Scx.,1 



, (29) 



G«a(0,A)| n=3 - fi _ 5U/2 _ x ■ /i _ 5C //2_2A 

°~a,,3 ^a.4 



H-5U/2-3X //-7Z7/2-4A 



, (30) 



In connection with the above observation, we note that 
by choosing the perturbation Hamiltonian in such a way 
that the state |<f„. So ), with s E {1, 2, . . . , Q }\{1}, 

would be the GS of T-Lx, we would have obtained expres- 
sions for lim>4o lirng-^oo Sf acc (0, A) (note the order of the 
limits), corresponding to different values of n, similar to 
those given in Eqs. (28) - (31), except for a non-trivial 
permutation of the 5 a ,j, j = 1> 2, 3, 4, in the numerators, 
which can be readily inferred from the expressions for 
the matrices {A a (n)||n} given in Eqs. (23) - (25) above. 
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This is relevant, in that it makes explicit that the zero- 
temperature limit of the thermal single-particle Green 
function ^(w), corresponding to H and [i £ (/i~,/i+), 

amounts to the arithmetic mean of {G^,(w)||s}, (^) 
single-particle zero-temperature Green functions, each 
corresponding to an n-particle GS of *H in {|#n,s)||s} 
singled out through introducing an appropriate pertur- 
bation Hamiltonian and taking the limit of the coupling- 
constant A of this perturbation approaching zero sub- 
sequent to effecting the zero-temperature limit (see the 
remark following Eq. (C.23) in Rcf. [2]). Thus, from the 
expressions in Eqs. (28) - (31) one trivially obtains that 



lira ^ a(0 ,0)| n=1 = -{ 

[3— >oo 4 l 



l r l 



(35) 



H - U/2 [i - 3U/2 . 

& ^« (0 < 0)l - = \{j^m + v^ml (36) 




} 

i 

/3->oo 



lim SP««(0,0)| n=3 = i{ 



\i - 3*7/2 /x - 5C//2 
3 1 , 



,(37) 



H-5U/2 H-7U/2P 
= \ u/2 > (38) 



where the integers 1 and 3 in the numerators are in fact 
multipliers of (<5«,i + 5 a ,2 + <5«,3 + ^a,i) = 1, the latter 
identity applying on account of a £ {1, 2, 3, 4} in the case 
of N = 4. 

By identifying the [i in the expressions in Eqs. (35) - 
(38) with Un, one readily obtains that 



lim lim Sf aa 



(o,A)L =t/ , 



U 



£{1,2,3,4}, (39) 



in full conformity with the expression in Eq. (13) of 
Ref. [1], which is specific to fi = Un, specialized to the 
case of N — 4. For the case of TV = 4, we have thus 
explicitly demonstrated that 

lim lim^aoiiui , A) # lim lim A), < n < N. 

fl^oo AJ.0 A J.0 ,3^oo 

(40) 

One can easily convince oneself that this result is not 
specific to N — 4, but applies for all N. The fact 
that despite this result the LT proves to apply for 
lim ( g_ ) . 00 lim^o A) when N is even, n = N/2 

and /it = (fi~ + /.t+)/2 = Un, is a direct conse- 
quence of particle-hole symmetry, whereby the latter 
zero-temperature single-particle Green function is equal 
to zero at u = (cf. Eq. (39), as well as Eqs. (9) and 
(C242) in Ref. [15]). It is important to realize that each of 
the zero-temperature Green functions contributing to the 
function lim^oo ^ aoc (0, 0) in Eqs. (35) - (38) satisfies the 

(s) 

requirement of the LT. In other words, with N^ denoting 
the LN corresponding to G^, (uj), the zero-temperature 
single-particle Green function corresponding to the n- 
particle GS |^ n ,s)j the thermal average of {N^ ||s}, in 
the zero-temperature limit, satisfies the LT. 



III. SUMMARY AND DISCUSSIONS 

Above we have first indicated that in the cases where 
the n-particle GS of a system is not unique, the zero- 
temperature limit of the thermal single- particle Green 
function Sf aa '(cj) corresponding to a fx € has no 

place in the formulation of the LT [2] , it being an aver- 
age (explicitly, arithmetic mean) of the zero-temperature 

Green functions {G^,(w)||s} corresponding to the set 
of n-particlc GSs {|^„ jS )||s}. The failure of the LT for 
the SU(iV) model under discussion for a general n, first 
established by DPK [1], is a direct manifestation of this 
fact. For the specific case of N = 4, we have explicitly 
shown how removal of the degeneracy of the GS energy 
by an infinitesimal amount leads to full restoration the 
LT. The transparency of the underlying calculations lead 
one to conclude that this result is general, applying to 
arbitrary values of N. The validity of the LT in the case 
of n = N/2, for TV even, and the chemical potential (i 
equated with the zero-temperature limit of (assuming 
that 0(0) = 1/2), despite the degeneracy of the relevant 
GS energy, is accidental, a consequence of particle-hole 
symmetry. 

From the experimental perspective, the GS energy of 
a system need not be degenerate in order for the LT to 
be observed as violated, as a consequence of the thermal 
energy T not being sufficiently small with respect to the 
lowest-lying single-particle excitation energies; for insuf- 
ficiently small T , the thermal single-particle Green func- 
tion Sf aa '(w) will necessarily contain considerable con- 
tributions corresponding to excited states, which, de- 
pending on the system under investigation, can prove 
destructive to the LT. For N = 4 and a fixed value of 
A > 0, we are able numerically to follow the development 
of Sf aa (w,A), Va. £ {1,2,3,4}, as a function of (3 and 
show how the LN N h as expressed in terms of ^ aa (0, A), 
corresponding to cither a fixed value of fi £ if-n'f-n) 
or u. = fip, deviates from n for /3 less than a critical 
value /3 C (A), and how N h coincides exactly with n for all 
f3 > /3 C (A). We provide the Mathematica® code underly- 
ing these calculations through an ancillary notebook file 
associated with the present text. 

We conclude this summary by stating that to our best 
knowledge to this date no case has come to light indica- 
tive of the failure of the Luttinger theorem under the 
conditions for which it has been deduced. 



Appendix A: Comments on some observations by 
Dave, Phillips and Kane [1] 

Some comments of general interest on some of the ob- 
servations by DPK in Ref. [1] are in place. For conve- 
nience of reference, we enumerate these comments. 

(i) The expression in Eq. (22) of Rcf. [1], with Ji and 
I2 as defined herein, is an identity, as emphasized in §4 
of Rcf. [2] (sec in particular Eq. (4.13) as well as §6.2.3 
of this reference). Consequently, reproducing the result 
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I± + I2 = n solely signifies correctness of the underlying 
calculations, and nothing more. 

(it) The singularity of the self-energy S aa >(uj) along 
the real 10 axis as deduced by DPK [1] (from the zero- 
temperature limit of the thermal single- particle Green 
function Sf aa '(u;) for [i € //+)), is a direct conse- 
quence of the non-uniqueness of the n-particle GS of the 
SU(N) model under consideration for < n < N. Be- 
fore elaborating on this statement, it is interesting to 
note that the self-energy E(oS) in Eq. (17) of Ref. [1] 
is identically vanishing for n — 0, N, the values of n for 
which the GS energy is non-degenerate. Here £{oj) is the 
short- hand notation for the diagonal element I7 aa (a;) of 
the self-energy, which proves to be independent of a, in 
conformity with the expressions on the right-hand sides 
of Eqs. (35) - (38) in the main text of the present Com- 
ment. 

Comparing the expressions in Eqs. (28) - (31) with 
their counterparts in Eqs. (35) - (38), one observes that 
the Green functions in the latter expressions, in contrast 
to those in the former ones, are independent of the index 
a, a fact fully attributable to the non-uniqueness of the 
n-particle GSs of the SU(iV) for < n < N. One readily 
verifies that the zero-temperature self-energy £ a0L {uj, A), 
Voc G {1,2, ... ,N}, associated with G aa (w,A), A > 0, 
and < n < N depends not only non-trivially on a, 
but also is independent of u> [16]. It is therefore analytic 
everywhere on the complex u plane. The independence 
from oc of the E aol (uj) as determined by DPK [1] and the 
simple pole of this function on the real co axis are inter- 
dependent and follow directly from the non-uniqueness 
of the n-particle GS of the SU(iV) for < n < N. 

For completeness, the behaviour of the self-energy 
along the real 10 axis plays no vital role in the proof of the 
LT. For this theorem it is however vital that So.cc' ( w ) be 
analytic anywhere away from the real axis of the complex 
ui plane [17], a fact emphasized repeatedly in Ref. [2] (see 
§§2.1.2, 5, 6.2 and Appendix B herein; see in particular 
the discussions related to the expression in Eq. (B.54)). 
The only point on the real energy axis where the be- 
haviour of £ 0,0.1 (u>) is to be given attention to is u = 
(according to the convention of Ref. [2] [§3.0.1], u) = fi, or, 
more precisely, z = fx). Considerations in Ref. [2] make 
explicit that insofar as the LT is concerned, from w = no 
problem arises (see §5.3.13, in particular the discussions 
centred on the expression in Eq. (5.61) of Ref. [2], as well 
as Appendix D herein). Interestingly, for < n < N, the 
singularity of the £(ui) as calculated in Ref. [1] is a simple 
pole located at lj — eo — M; where eo is defined in Eq. (18) 
of Ref. [1]. On identifying fi with = +/z+)/2, for 
U > this pole is located at u — only if N is even and 
11 = N/2. Remarkably, for this particular case the LT in 
terms of the N h presented in Eq. (5) is satisfied. 

(Hi) With reference to the data displayed in Fig. 2 of 
Ref. [1] concerning the "apparent doping x FS inferred from 
the Luttinger surface reconstruction as a function of the 
nominal doping x in LSCO and Bi-2212" [1], we are not 
in a position to make a definitive statement about these 



in the absence of any information regarding the method 
by which x FS versus x has been inferred from the exper- 
imental data of He et al. [18] and Yang et al. [19]. For 
completeness, insofar as the LT is concerned, in these 
references only very brief remarks are to be found, which 
furthermore are suggestive of qualitative agreement with 
the LT: on p. 7 of Ref. [18] one reads: "the k F separation 
shows a systematic increase with doping, which is con- 
sistent with the increase of Luttinger's volume of the FS 
[Fermi surface] and the corresponding shift of the node 
position away from fc AF " (A; AF is defined in the caption of 
Fig. 1 of Ref. [18]), and in Ref. [19]: "One does not need 
to invoke discontinuous Fermi "arc"s to describe the FS 
of underdoped Bi2212 and Luttinger's sum rule, properly 
understood, is seen to still approximately stand". 

In spite of the above observations, two remarks are in 
place. Firstly, experimental determination of the 'vol- 
ume' of the Fermi sea is complicated by an aspect that 
to our best knowledge was first spelled out by Essler and 
Tsvelik [20], and Konik, Rice and Tsvclik [21]. In Ref. [2, 
§2] we have highlighted the problem and remarked that 
the Fermi sea of a metallic GS not being necessarily a 
closed set of k points, knowledge of a Fermi surface is in 
general not sufficient for determining the corresponding 
Fermi sea (see the discussions following Eq. (2.7) in §2.1 
of Ref. [2]). 

Secondly, in §6.4 of Ref. [2] we have exposed the mech- 
anism by which in the case of strongly-correlated met- 
als (for which the single-particle momentum distribution 
function n CT (fc), pertaining to particles with spin index c, 
strongly deviates from a unit-step function that is char- 
acteristic of the GSs of non- interacting fermions), some 
regions of the k space that are in reality external to the 
underlying Fermi sea, can be mistakenly identified as re- 
gions internal to this sea. As discussed in detail in Ref. [2, 
§6.4], this mechanism underlies the erroneous observa- 
tion by Grober, Eder and Hanke [22] regarding the LT 
(see in particular Fig. 12 in Ref. [22] and note that this 
figure displays the "measured" volume of the Fermi sea 
versus the hole concentration 1 — n, defined in relation to 
the half- filled, n = 1, state of the single-band Hubbard 
Hamiltonian - here, in two dimensions). 

(iv) With reference to an earlier work by Stancscu and 
Phillips [23], to which DPK [1] refer, the calculations in 
this work have been performed under a number of ap- 
proximations (e.g. 'the two-site approximation'), spelled 
out in §11 of Ref. [23] (see specifically pages 5 and 6 
herein). Aside from this, the computations reported in 
Ref. [23] correspond to the w-dependent functions that 
have been "discretized on a grid of N = 8192 points from 
^min = — 20t to uj max = 20t", with "all the convolutions 
involved in the calculation of self-energies" having been 
evaluated through "using a fast Fourier transform algo- 
rithm." [23]. Here t denotes the nearest-neighbour hop- 
ping integral, and N should not be confused with the N 
of the SU(7V) model. Further, "The procedure converges 
for temperatures above T = 0.02t at finite doping and 
T = 0.08t at half-filling, although convergence problems 
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occurred below T = O.li." [23, p. 9]. The Fermi-surfacc 
geometries displayed in Fig. 18 of Ref. [23] correspond to 
band fillings n = 0.3, 0.668, 0.791 and 0.97, and U = 0, 
8t and lOOOi. 

In Sec. II B 2 of the present Comment we have dis- 
cussed the adverse effect of non-zero temperatures for 
the validity of the LT. In the following we shall therefore 
focus on other aspects that unquestionably render the 
observations in Ref. [23] regarding the LT unreliable. 

In §6.3 of Ref. [2] we have extensively discussed the 
detrimental consequences of using insufficiently large cut- 
off energies / frequencies (denoted by — E and E in 
Ref. [2]) for the self-energy and in particular the LT. In 
doing so, we have concentrated on some relevant works 
by Schmalian et al. [24, 25]. In the calculations reported 
in these references, E = 30i, where t = 0.25eV, and 
N = 4096. Clearly, the cut-off frequencies w m j n and w max 
in the calculations by Stanescu and Phillips [23] are of the 
same order of magnitude as respectively —E and E. The 
lowest temperature considered in Refs. [24, 25] amounts 
to T = 63 K, which, making use of the relationship 1 
cV « 1.16 x 10 4 K, amounts to approximately 0.022£. 
The details underlying the numerical results presented in 
Ref. [23] are not as exhaustively described as those un- 
derlying the numerical results presented in Refs. [24, 25]. 
It is therefore not possible for us to be as specific with 
regard to the former numerical results as we have been 
with regard to the latter ones in Ref. [2, §6.3]. However, 
from the data for the density of the single-particle states 
(DOS) as presented in Fig. 11 of Ref. [23], one can in- 
fer that a combination of insufficiently small value of T 
and insufficiently large value of |w m i n | = w max plays a 
significant role in the violation of the LT as observed in 
Ref. [23] , with the latter condition being the more impor- 
tant of the two. Sec Fig. 2 of Ref. [24], to be compared 
with Fig. 4 of Ref. [2] (p. 107), and consider the details 
centred around Eqs. (B.42) - (B.46) in Appendix B of 
Ref. [2]. 

We note in passing that, violation of the exact prop- 
erty lm[E(k, 0)] = 0, Vfe, is the fundamental reason un- 
derlying the observation by Maier, Pruschke and Jar- 
rell [26] of the breakdown of the LT in the case of 
the single-band Hubbard Hamiltonian in two dimensions 
and at low concentration of holes introduced into half- 
filled GSs [2, §6.3.6]. For the relationship between the 
DOS, X^fceiBZ A(k\ u), where A(k; ui) denotes the single- 
particle spectral function, and Im[U (k,uj)], the reader is 
referred to Eq. (B.43) in Ref. [2]. 

The case of U = lOOOt (in general, large values of 
U/t) as considered by Stanescu and Phillips [23] deserves 
special attention, this in the light of the considerations 
in Appendix D of Ref. [27]. With w m j n = -20f and 
^max = 20i, the choice of U = 1000£ takes one artifi- 
cially to the region of Extremely Correlated (EC) limit 
[27, Eq. (D.2)], whereby breakdown of the LT becomes 
a certainty, on account of some non-vanishing boundary 
terms that appropriately have not been taken account 
of in the derivation of the LT (since in this derivation 



E = oo). In the light of the explicit use in Ref. [27, 
Appendix D] of the single-particle Green function corre- 
sponding to the atomic limit of the Hubbard Hamiltonian 
[6, §5] , we draw the attention of the reader to the remarks 
following Eq. (10) in the main text of the present Com- 
ment. 

To summarise, the breakdowns of the LT as observed 
in Ref. [23] are to our best judgement numerical arti- 
facts. On account of the detailed considerations in §6.3 
of Ref. [2], we are fully confident that on redoing the 
calculations with sufficiently large values of — Wmi n and 
w max (large in comparison with U), the LT will prove to 
be valid. 

(v) In contrast to the statement by DPK [1, p. 4], in 
Ref. [28] , no "systematic deviation" from the requirement 
of the LT has been observed for the Hubbard model. 
In the concluding section of Ref. [28], p. 5, one reads: 
"For the Hubbard model on the ID chain and 2D square 
lattice, both for t' = and t! ^ 0, we do not find a 
clear-cut violation [type (iv)] of the LSR [Luttinger sum 
rule]." 

As regards the violation of LT for the t-J model, for 
which this violation has indeed been observed in Ref. [28] , 
there is no a priori reason why the LT should be valid 
for this model. In fact, the theoretical considerations in 
Ref. [27] make explicit that this violation is inherent to 
the t-J model. In this connection, we point out that the 
criterion employed in Ref. [29] (while emphasizing its lim- 
itation) for determining the Fermi surface of the metallic 
states of the t-J model, namely that the Fermi surface 
corresponding to fcrmions of spin index o~ were the locus 
of the k points on which n CT (fc) = |, where n a (k) denotes 
the underlying GS momentum distribution function, is 
unfounded [30] (see in particular the Appendix herein). 
The observation in Ref. [29] of the violation of the LT, 
based on the latter defining equation for the Fermi sur- 
face, can therefore not be viewed as reliable and thus 
conclusive. 

For completeness, as indicated earlier by a number or 
authors (such as Chao, Spalek and Oles [31], Eskes et 
al. [32], Dagotto [33], and Eskes and Eder [34]; see also 
Ref. [7, Ch. 5]), the strong-coupling limit of the Hub- 
bard Hamiltonian involves in addition to the t-J Hamil- 
tonian some three-site terms, proportional to t 2 /U. In 
identifying the t-J Hamiltonian as the strong-coupling 
limit of the Hubbard Hamiltonian, one therefore dis- 
cards the latter terms, a practice that has no rigorous 
justification away from half-filling. In this connection, 
we should emphasize that even at half-filling one re- 
tains an active three-site term, arising from the anti- 
commutation [iS, H®) (in the notation of Ref. [7]) asso- 
ciated with the canonical transformation through which 
the strong-coupling effective Hamiltonian H c ff is deduced 
from the Hubbard Hamiltonian H [7, p. 210]. Note that 
in contrast to the Hubbard model for fermions for which 
the Fock space of each site is spanned by four states, |0), 
|t)j II) an d | tl)j f° r the t-J model the Fock space of 
each site is spanned by only three states, with the doubly- 
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occupied state | fi) excluded. 



Appendix B: Discussion of a recent work by Sakai et 
al. [35] 

In a recent publication regarding calculation of the self- 
energy S of the Hubbard Hamiltonian on finite clusters 
in two space dimensions, with the aid of the 'cellular 
dynamical mean- field theory', Sakai et al. [35] have ob- 
served an indication "strongly" suggestive of the failure of 
the LT. On the basis of the observed "cluster-size depen- 
dence of £ remaining around (tt,0) and (tt, 7r)", making 
it "difficult to make definitive statements", the authors 
have postponed their definite verdict on the violation or 
otherwise of the LT until after larger clusters will have 
been studied [35, §111. C]. 

A noteworthy aspect of the computational results by 
Sakai et al. [35] is that their calculated self-energy is 
£(k,iu>o), where ujq = tt//3, the Matsubara frequency 
Lu m = (2m + 1)tt//3 corresponding to m = 0; the temper- 
ature T = 1/ j3 in the calculations of Ref. [35] is chosen to 
be equal to 0.06t, where t denotes the nearest-neighbour 
hopping integral (other parameters in these calculations 
are t' = —0.2t and U = 8t, where t' is the next-nearest- 
neighbour hopping integral, and U the on-site interaction 
energy). Partly, however not essentially as we shall dis- 
cuss below, result of cjo riot being infinitesimally 
small, the imaginary part of £(k,iujo), as depicted in 
the right-most panels of Figs. 6 and 7 in Ref. [35], is 
non-vanishing (in fact, it is unusually large - see the fol- 
lowing paragraph). This is inadmissablc from the per- 
spective of the LT, since for the exact self-energy in the 
zero-temperature limit one has lm[U(k, 0)] = 0, Vfc (see 
§2.1.2 and Appendix B in Ref. [2]; compare also with the 
numerical results for the real and imaginary parts of the 
self-energy in Fig. 1 of Ref. [36] and note that these re- 
sults correspond tot = 1 and U = 2). See the discussions 
under point (iv) in appendix A. 

On general grounds, one expects the variation of 
Im[S(k, iwo)] as function of k to be on the scale of loq 
(cf. Eqs. (B.55) - (B.64) in Ref. [2]), which in the case 
at hand is approximately equal to 0.19i, in contrast to 
Re[27(fe, iwo)] whose variation is expected to be on the 
scale of, say, U. The numerical results presented in 
Figs. 6 and 7 of Ref. [35] show that for k varying over the 
underling fe-spacc, both Kc[S(k, iu>o)] and Im[S(k,iuio)] 
vary over intervals whose widths are of the order of U. 
This behaviour is inexplicable to us. Hence, it is essen- 
tially ruled out that use of S(k,iujo), instead of S(k, 0), 
could be the reason for the numerical results in Ref. [36] 
strongly suggesting violation of the LT. The employed 
computer codes need to be inspected for error (s). 



Appendix C: A Reply to two Replies by Kokalj and 
Prelovsek [37] 

In this appendix we respond to two Replies by Kokalj 
and Prelovsek (hereafter KP) [37] on the Comments in 
Ref. [15] (co-authored by A. M. Tsvelik) and Ref. [38]. 
Below, the numbers 1), 2), . . . refer to the itemized replies 
in Ref. [37] numbered thus. 



1. Concerning the Comment by Farid and Tsvelik 

[15] 

1) The remark in Ref. [15] "Neither Stanescu, Phillips 
and Choy [11] nor KP [12] appear to have appreci- 
ated this fundamental aspect.", to which KP object in 
Ref. [37], has direct bearing on the following statement 
made in Ref. [39, p. 1]: "On the other hand, there are 
several indications that LSR might be violated within 
the MH [Mott-Hubbard] insulators in general. 8 " 10 " 

2) - 

3) As for the supposed "claim" made in Ref. [15], in 
Ref. [39] KP state: "According to the LSR, for the spin- 
less model at n = 1/2 one should generally have k L = tt/2 
if the topology of the electronic band is not changed qual- 
itatively (which could happen, e.g., for t' > 0.5i). From 
Fig. 2 we note that k L is indeed near tt/2; however, even 
without finite-size scaling [our emphasis] a small devia- 
tion k h 7^ tt/2 may be observed for V > 4i." 

In Ref. [15] we have made the following two main ob- 
servations: 

(i) For N = 26, where N denotes the number of lattice 
sites, 7r/2 does not belong to the underlying k space. For 
this very reason, for N — 26 the fc L corresponding to 
half-filling cannot have coincided with tt/2. In this light, 
and "without finite-size scaling", fc L ^ tt/2 cannot have 
signalled violation of the LT, contradicting the explicit 
claim to the contrary by KP [37, 39] (see above) . 

Although indicated earlier [15], we emphasize that the 
LT (or the Luttinger sum rule, LSR) states that 7V L = N c 
[15, Eq. (2)], where N h , the LN, is defined in Eq. (1) of 
Ref. [15], and N e denotes the number of particles in the 
GS under consideration. From Fig. 2 of Ref. [39] one 
clearly observes that for all values of V/t shown (barring 
for the moment the value of V/t = 8 on which KP focus 
in their pertinent Reply [37]), at exactly 13 fc-points, out 
of the total of 26 available Appoints, one has G(k, /i) > 0, 
in perfect agreement with the LT, given the fact that the 
GSs under consideration are half-filled. 

We remark that from Fig. 2 of Ref. [39] it appears 
that the solid line (drawn in red) , presented as depicting 
the locus of k L as a function of V/t, is determined by 
taking the zero of the straight line connecting the val- 
ues of G(k, 0), corresponding to a given value of V/t, at 
k = 67r/13 and k = 7tt/13. Since the centre of the solid 
dot representing the value of the G(k,0) corresponding 
to V/t = 8 and k = 7tt/13, denoted by ki in Ref. [37], 
is visibly to the right of the middle of the red solid line 
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depicting k L , we conclude that the "G(ki,ui — 0) ~ 0" of 
Ref. [37] in fact represents a small however negative value 
for G(ki,tJ = 0), conform the requirement of the LT. Be 
it as it may, we feel frustrated at the fact that KP in 
their pertinent Reply [37] do not disclose the actual nu- 
merical value of G(ki,uj = 0) and only suffice to present 
the relationship "G(ki,LJ = 0) ~ 0". It should be noted 
however that in the absence of any published details re- 
garding the accuracy of the numerical results presented 
in Ref. [39] , even a positive but sufficiently small value of 
G(ki , ui = 0) cannot be taken as signifying failure of the 
LT. 

As for the numerical results corresponding to A = 30, 
to which KP refer in Ref. [37], since these are not pre- 
sented in Ref. [39] (nor are they in public domain), our 
Comment [15] could by no means have direct bearing on 
them. Further, as emphasized in Ref. [2, §2.5], for finite 
systems the k points at which G(fc,0) ~ have to be 
treated with care in taking the zero-temperature limit of 
the expression for the LN N h . 

(ii) In the thermodynamic limit, the GS of the half- 
filled one-dimensional t-t'-V model for small values of 
t'/t and V > V c (t,t'), where V c (t,t') ~ 2t for t'/t 0, 
is a charge-density-wave (CDW) state and the gap in 
its single-particle excitation spectrum is not a MH gap, 
but an ordinary gap arising from translational symmetry 
breaking [15]. In the light of this observation, in §IV of 
Ref. [15] we have explicitly shown that in the insulating 
phase of the model under consideration the LT trivially 
holds. 

4) The toy model employed in Ref. [15] is not essential 
for the observations in this work. It is merely intended to 
show that for insufficiently large values of A, finite-size 
scaling can yield very inaccurate values for fc L . 

5) The argument presented by KP in Ref. [37] only 
further discredits their use of finite-size scaling method 
for arriving at a conclusion with regard to the validity or 
otherwise of the LT for the t-t'-V model in the thermo- 
dynamic limit. 

Regarding the statement in Ref. [37], that "Also it 
should be pointed out that the long-range order does 
not necessarily invalidate the LSR but merely needs a 
redefinition of the latter as correctly pointed out in the 
Comment", in §IV of Ref. [15] we have explicitly dealt 
with this very redefinition, due to Luttinger [3, Eq. (95)]. 



2. Concerning the Comment by Farid [38] 

The statement by KP in Ref. [37] that we merely on 
"p. 3 (after Eq. (28))" of Ref. [38] had 'claimed' "fallacy 
in the reasoning" of Ref. [40] , suggests that we must not 
have been sufficiently clear and exhaustive in our argu- 
ments in Ref. [38] . Below we attempt to compensate for 
these shortcomings and in doing so present some new in- 
sights that we have gained since we published Ref. [38] 
in September 2009. 

1) We do not make any assumption of the kind sug- 



gested by KP [37] in the analysis of Ref. [38]. In this 
reference we merely argue that cannot be justifiably 
identified with fc L , thereby to arrive at a definitive con- 
clusion with regard to the validity or otherwise of the LT. 
It is relevant to recall that the vector k^"\ with v 6 N, 
was in explicit form first introduced in Ref. [38] , and not 
in Ref. [40]. 

2) We do not object to the expansion in powers of 
fx = fj, — U/2, to be distinguished from the expansion 
in powers of t/U (see later). As a matter of fact, just 
immediately before Eq. (5) of Ref. [38] we directly refer to 
the expansion in powers of jl herein as "exact". Further, 
following Eq. (6) of Ref. [38] we fully specify the interval 
of jl over which the expansion in Eq. (5) is valid, and 
nowhere in Ref. [38] do we indicate or even imply that ji 
may be outside the latter interval; the chemical potential 
\x corresponding to the A-particle GS of the system under 
consideration does not allow for this to be the case. 

3) The problem with the considerations of Ref. [40] is 
two-fold, fundamental and technical. Before describing 
these, below for convenience we consider ty = t 
for i and j nearest neighbours, and = otherwise. 
In the cases where the latter conditions are not met, the 
considerations are to be based on the quantity A, instead 
of t, where A is introduced through the identity Uj = 
A nj , in which J2j M 2 - 1 [41, Eq. (2.4)]. 

Fundamentally, insofar as the LT is concerned, t/U = 
is a priori not an appropriate limit around which to per- 
form a finite-order expansion in powers of t/U . This is 
rooted in the fact that the GS of the Hubbard Hamilto- 
nian in the atomic limit [6, §5] is not unique (Sec. II A, fol- 
lowing Eq. (10)). Consequently, the LN N h as calculated 
in terms of the zero-temperature limit of the thermal 
Green function ^o-o-'(O), corresponding to fx € 
is in general not equal to n. From the perspective of the 
LT, t/U = is thus the singular point of the strong- 
coupling expansion of the single-particle Green function. 
Here cr, & 6 {til} denote spin indices for spin-i parti- 
cles. 

Technically, there are two distinct aspects to be con- 
sidered. Firstly, close inspection of the work by KP in 
Ref. [40] reveals that the underlying calculations are ex- 
actly at the t-J level of the strong-coupling expansion of 
the Hubbard Hamiltonian [7, Ch. 5]. It is therefore not 
surprising that a strict adherence to the adopted low- 
order formalism of Ref. [40] should signal breakdown of 
the LT (see entry (v) in appendix A as well as Ref. [27]). 
For clarity, on disregarding the rest terms in the expres- 
sions for the functions Mq (fe) and Aff(fe) in Eq. (11) 
of Ref. [39], one is left with functions that arc spe- 
cific to the t-J model at half-filling, or the Heisenberg 
model (up to an unimportant additive constant due to 
the density-density interaction term in the t-J Hamilto- 
nian) [7, §5.1.5] . 

Secondly, assuming that the results G^(k { L 1] ; U/2) = 
0(t/U 4 ) and G(bO)(fe L ; M ) = Q(t/U 2 ), deduced by KP 
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[37], are exact (see Sec. C2a), whereby 

G(°- 2 >(fc<V/2) / 1 x 
G( 1 '°)(fc L ;/i) It/V' 



(CI) 



the asymptotic expression in Eq. (16) of Rcf. [38] for k^ 
in terms of k L becomes inapplicable (for the considera- 
tions at hand, v = 1). This follows from the fact that by 
the same reasoning that G^ 1,0 ) (fc L ; fx) = 0(t/U 2 ), one has 
Gf(m,o)( fct . n) = 0{t/U 2 ) [38, Eq. (4)] for any finite value 

of to, Sec. C2a. As a result, so long as \k^ — k L \ is not 
sufficiently small, at the largest of the order of (t/U) 2 (we 
assume the lattice constant a = 1), the sum with respect 
to to in Eq. (15) of Rcf. [38] cannot be approximated by 
its summand corresponding to to = 1 . 

With G^ 2 \k { h 1} ;U/2) = 0(t/U 4 ) [37], Sec. C2a, one 
is thus led to the conclusion that, on account of Eq. (12) 
in Ref. [38], to leading order in t/U, as t/U — > 0, the 

equation to be solved for fc L is G(fe£ ,fi) = 0, viewed 
as an implicit function of fe L . Expressed differently, fe L 



is the solution of the equation G(k\fx) = under the 
subsidiary condition that to order t/U 2 (inclusive) the 
vector k^ also satisfies that latter equation. To clar- 
ify, by definition fe^ = &l ; n is the exact solution of 
G v (k;fx) = [38, Eq. (13)]. The asymptotic expression 
for fc^ in Eq. (16) of Ref. [38] being deduced by equat- 
ing the asymptotic expression for G(k^ ; fx) on the right- 
hand side of Eq. (14) of Rcf. [38] with the aforementioned 
summand corresponding to to = 1 in the expansion of 
G(fe£ ; fi) in powers of (fc£ — k L ), it neglects an infinite 
number of terms, each of which through {G^ m ' ' (k L ; fx)} 
scales to leading order like t/U 2 , Sec. C2a, in favour 

of a term, i.e. G^' v+1 ^ [k^ ; 17/2), that for t/U to 
leading order scales like t/U" +3 when v is odd and like 
\/U u+2 when v is even, Sec. C2a. 

Following the above considerations, for a given to 
leading order in t/U the Luttinger vector k L is the solu- 
tion of the following asymptotic equation (cf. Eq. (15) in 
Rcf. [38]): 



£ G(m ' 0) J*<>ri (fcW - k h ) m = G(feW; M ) - G(fc L ; M ) ~ 0. (C2) 

m— 1 

Clearly, this equation is trivially satisfied for k h = fc£ , the solution on which KP's [37, 40] conclusion regarding the 

LT is based. Dividing both sides of Eq. (C2) by (&l — k L ), and assuming that G^ (k h ; fx) ^ for to = 1,2, one 
can express the resulting equation in the following appealing form: 



GWifan) r , G ('"+ 1 '°)(fc L ; / ,) (A>> - fc L ) m 



G(2.o)(fe L ; M )l ^ 2 Gd-o)^;//) (to + 1)! J' 



The right-hand side of this asymptotic expression de- 
pending of fc L , A; L is to be determined iteratively. Only for 
sufficiently small \k^ — k L \ may one identify the k L and fe L 

on the right-hand side by respectively and k^ , and 
further employ the following approximate asymptotic ex- 
pression: 



, G^\k^ ]f x) 
G( 2 >°)(fc^;/x) 



(C4) 



Note that by the considerations of Sec. C 2 a, one has (for 
the lattice constant a = 1) G (1 ' 0) (fc L ; /i)/G (2 '°) (fc L ; /x) = 
0(1) and G( m+1 '°)(fcL;M)/G (1 ' 0) (fcL;M) = 0(1), Vto > 1. 
In Fig. 1 we display a contour plot of the function 

on the left-hand side of Eq. (C2) divided by (kl — k L ). 
The function G(fc; /x) employed in the calculations is that 
given in Eq. (C9) for to = 0, which is exact to order t/U 2 
in the region t/U — > 0, provided that the underlying 
functions M^(k) and Mf(fc) arc calculated exactly to 
leading order in t/U. The functions M^(k) and M^(k) 
employed in the present calculations are those presented 



in Eq. (11) Ref. [40], which, as we have indicated earlier 
in this section, are specific to the t-J model at half-filling, 
resulting in the expression for G(k; fx) as given in Eq. (16) 
of Ref. [40] (Eq. (20) in Rcf. [38]). Following this obser- 
vation, it should not come as a surprise that within the 
present approximate framework, Eq. (C2) has no other 

solution than k h = fc£ ■ To appreciate the significance of 
the shortcoming as arising from the t-J approximation of 
the Hubbard Hamiltonian in the strong-coupling region, 
one should consider the contour lines in Fig. 1, showing 
that a rigid upward shift of the underlying function by 
merely approximately 0.006 is sufficient for producing a 
k h that by the required amount (from the perspective of 
the LT) is larger than k[ . See in particular the con- 
cluding paragraph of Sec. C 2 a. For orientation, with 
ft = t(0.1 + 6.8t/U), for n in the positive x direction of 
the underlying 1BZ (see Eq. (25) in Ref. [38]), one ob- 
tains k^/ir « 0.725, to be contrasted with k h /n « 0.849, 
deduced by requiring satisfaction of the LT. To a good 
approximation, these two values clearly correspond to a 
point on the contour in Fig. 1 labeled —0.00568. 
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FIG. 1. Contour plot of the left-hand side of Eq. (C2) di- 
vided by (fct — fa.) on the k^-k L plane for both fcf^' and fc L 
along the positive direction of the horizontal axis of the 1BZ 
as specified in Eq. (25) of Ref. [38], indicative of k L = fc£ 
as being the only solution of Eq. (C2). As we discuss in the 
main text, this result was to be expected, on account of the 
function G(k; /i) employed here corresponding to the t-J ap- 
proximation of the strong-coupling expansion of the Hubbard 
Hamiltonian, Sec. C 2 a. The data displayed here correspond 
to t = 1 and U = 40, with the lattice constant a identified 
with unity. 



For completeness, and with reference to the expression 
in Eq. (16) of Ref. [38], we note that the order of mag- 
nitude that we have presented in Eq. (29) of Ref. [38] 
for the coefficient of the latter equation as specialised to 
the case of v = 1, is based on the Kramers-Kronig re- 
lation and the general property of Im[G(fc; uj ± i0 + )] for 
insulating GSs as u> passes through the band edges (see 
reference [16] in Ref. [38]). Further, since in Ref. [37] KP 
refer to Eq. (29) of Ref. [38] (their Eq. (2)) as ambigu- 
ous, we draw the attention of the reader to Ref. [42]. At 
any rate, we concede that the expression in Eq. (16) of 
Ref. [38] has no root in the strong-coupling expansion of 
the single-particle Green function G(k; /i), in that it does 
not take into account the details that lead to the equal- 
ities G^ 1 ' ) (fe L ; n) = 0(t/U 2 ) and G^^{k^ ] -U/2) = 
0(t/U 4 ). It is to be noted however that whereas the 
latter equality is exact, the exactness of the former one 
is not rigorously established. For this, see the remarks 
in the paragraph following that containing Eq. (C8) in 
Sec. C2a. 

a. Expressions for G (0 ' n) (fc; 17/2) and G {m ' 0) (fe; /i) 

In order to avoid notational confusion, unless we in- 
dicate otherwise, in this section we make use of the ex- 
plicit expressions in Ref. [38], where we have appropri- 
ately credited Ref. [40] where credit has been due. 

From the expressions in Eqs. (8) and (11) of Ref. [38] 
one has 



G(k^) = J2(-Hr(±) {A/ r( fc)-(-l)<+«M+(fc)}, (C5) 



n=0 z=o 



where M^(k) are defined in Eqs. (6) and (7) of Ref. [40]. Since Ji = fx — U/2, one has d n /dfi n = d n /dfi", Vn, and 
further fl = for \i = U/2. Thus, from Eq. (C5) one trivially obtains that (cf. Eq. (4) in Ref. [38]) 



G<°-»>(fe; U/2) = (-ir(!)" +1 £ *4^(§)V<~ (fc) " (- 1 ) i+ " M ; + ( fc )} 

i=0 



/2\»+ir /2\«+i 2n fcs -I, n = even, , 

(-!)"«!(-) |M -(fc) - (-l)"M+(fe)} = n\(-) x I for - -> 0, (C6) 

n = odd, 



where in the last expression we have employed the first 
two equalities in Eq. (8) of Ref. [40]. Identifying nfe s , 
the GS momentum distribution function correspond- 
ing to particles with spin index s £ {t)4-}i with that 
corresponding to the Heisenberg model [32, Eq. (31)] 
[40, Eq. (11)] [7, Eq. (5.45)], one arrives at the result 
2n ks - 1 = 0{t/U) [43], specifically for k = fc[ 1} . Hence 
the equality C?( 0,2) (kP;U/2) = 0{t/U 4 ) as presented in 



Ref. [37]. 

Making use of the relationship in Eq. (3) of Ref. [38] 
and defining, in analogy with the G^ m ' n ^(k; fj,) in Eq. (4) 
of Ref. [38], 

Mf {m \k) = gj^M^ik), m = 0, 1,2, . . . , (C7) 
from the expression in Eq. (C5) one deduces that 
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G (m ' 0) (k;p) 



v E(-«" £ (I) (|) K-L m) (fe) - (-d'm^w}. 

n=0 Z=n ^ ' 



(C8) 



Considering the expression for ft = — {7/2 as presented 

in Ref. [40], namely // (0.19 ± 0.1)t + 6.8£ 2 /?7 (cf. 

Eq. (27) in Ref. [38]), one observes that G( m '°)(fe L ; /i) is 
a non-trivial function of and i. 

Some technical remarks concerning the equality in 
Eq. (C8), for m > 1, are in place. This equality has 
been deduced by m times commuting d/dk with infinite 
sums with respect to n and I in the expression for G(fc; fx) 
(this expression is up to a trivial re-indexing of the sum- 
mation with respect to I, similar to that in Eq. (C5)). 
For general fc-dependent summands, this is admissible 
only if the sums with respect to n and / on the right- 
hand side of Eq. (C8) are uniformly convergent for k in 
a neighbourhood of the k of interest, such as k = k L (as- 
suming that these summands are continuous functions 
of k throughout the last-mentioned neighbourhood and 
that the relevant sums corresponding to G 1 -'™" 1-0 ' (fc; /1) 
are convergent). For the cases where the summands of 
the sums with respect to n and I corresponding to the 



I 

function G < - m_1,0 ^(fc; fx), m > 1, are analytic functions of 
k in a neighbourhood of the k of interest, such as k = k L , 
the uniform convergence of these sums for k in the latter 
neighbourhood suffices for commuting d/dk with them, 
whereby G' 1 ™' - 1 (fe; fx) is expressible as in Eq. (C8) [44, 
§§4.7,5.3], [45, §46]. We note that the above-mentioned 
sums with respect to n being power series, they are 
convergent in the region specified in the text following 
Eq. (6) in Ref. [38]. In fact, in this region these sums 
converge uniformly with regard to the variable ft [44, 
§3.7], [45, §§50-52]. As regards continuity of G^ m '°\k; fx), 
m > 0, in its dependence on k, see Ref. [45, §53]. We note 
in passing that it is in general not permissible to differ- 
entiate asymptotic series [44, §8.31]. 

From the expression in Eq. (C8) one infers that the 
lowest-order contributions to G^ m '°'(fc;/i) for t/U — > 
arise from the terms corresponding to (n, I) equal to 
(0, 0), (0, 1) and (1, 1). For t/U —¥ one thus has 



G(-°)(fc; M ) ~ |{A/-^(fc)-A/ + ^(fc)} + (|) 2 {[Mr (m) (fc) + A/ 1 +(m) (fc)] - [M-^\k) + M^ m \k)]ft}. (09) 



Making use of the expressions in Eq. (11) of Ref. [40], one 
arrives at the result G (1 ' 0) (fc L ;^) = 0(t/U 2 ), presented 
in Ref. [37]. More generally, one has G (m < 0) (fc L ; fx) = 
0(t/U 2 ) for any finite value of m, including m = 0. 
Clearly, the validity of this result for m > 1 is depen- 
dent on that of the expression in Eq. (C8). For this, see 
the paragraph following the latter equation. 

We point out that since M "(fc) + M^(k) = 1 for all 
k, it follows that the function on the right-hand side of 
Eq. (C9) multiplying the ft is identically vanishing for all 
m > 1. Consequently, G^ m '°\k; fx) is to order 1/U 2 inde- 



pendent of ft for all m > 1. This is a manifest shortcom- 
ing of the functions {G^ m '°^(fe L ; fx)} calculated to leading 
order in t/U for use in the calculation of k h according 
to the expression in Eq. (C2), or that in Eq. (03). In 
particular with reference to the latter equation, we note 
that in general for / ~ a^e 2 + bfe 3 and g 



a 



as e — > 0, one has f/g ~ aj/a g + (bf/a g — afb g /a 2 )e 
as e — > 0. For comparison, owing to the exact identity 
M -(fc) + M+(fc) ee 1, G(°'°)(fe; n) ee G(fe; n) is to order 
1/U 2 linearly dependent on ft (cf. Eq. (16) in Ref. [40] 
and Eq. (20) in Ref. [38]). Consequently, ki , in partic- 
ular for v = 1, sensitively depends on the value of ft. 
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